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ABSTRACT 


The  application  of  numerical  optimization  techniques  to 
the  design  of  drive  shafts  is  demonstrated.  The  analysis 
investigated  the  effects  of  a  small  mass  imbalance  in  con¬ 
junction  with  the  rotation  of  a  shaft  with  synchronous 
whirl . 
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I.  INTRODUCTION 


Optimization  and  composite  materials  are  two  fields 
that  have  seen  considerable  developments  in  recent  years. 
Advancements  in  these  technologies  have  provided  the  design 
engineer  with  new  dimensions,  the  depths  of  which  are  yet 
to  be  explored. 

The  major  material  characteristics  that  are  advanta¬ 
geously  cultivated  in  composites  are  the  high  strength  and 
stiffness  to  weight  ratios.  The  weight  reduction  is  so 
significant  that  the  aircraft  and  aerospace  industries  have 
concentrated  a  major  portion  of  their  reserach  efforts  in 
this  field.  Investigations  in  improving  the  resistance  of 
composites  against  environmental  factors  such  as  temperatures, 
corrosion,  and  wear  are  continuing  and  findings  affirm  the 
contention  that  the  future  of  this  field  is  bright.  However, 
other  industries  are  slow  in  utilizing  these  materials  in 
basic  machinery  elements,  mainly  due  to  a  major  drawback — 
the  high  cost  of  manufacturing  and  processing.  It  is  in 
this  vein  that  the  importance  of  an  optimized  design  becomes 
apparent.  Advancements  in  fabrication  and  processing  of 
composite  materials  coupled  with  the  capability  to  optimize 
designs  may  well  lead  composites  into  a  more  competitive 
market. 
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With  the  ultimate  goal  of  optimizing  the  design  of  a 
drive  shaft  made  of  composite  materials,  a  research  program 
was  initiated.  This  report  is  the  first  in  this  investiga¬ 
tive  effort. 

Before  the  intricacies  introduced  by  the  directional 
material  properties  of  composites  are  included  in  the 
analysis,  it  was  deemed  necessary  to  probe  the  characteris¬ 
tics  of  the  design  of  an  isotropic  material.  This  will 
enable  one  to  foresee  the  trends  and  the  adaptive  measures 
necessary  when  the  composite  material  properties  are  con¬ 
sidered.  The  application  of  numerical  optimization  tech¬ 
niques  to  the  design  of  an  isotropic  drive  shaft  is  demon¬ 
strated  here.  The  effects  of  a  small  mass  imbalance  on  the 
bending  moment  and  deflection  of  a  shaft  with  synchronous 
whirl  are  investigated.  Mainly,  however,  the  objective  is 
to  develop  an  analysis  procedure  that  may  be  adapted  to  the 
peculiarities  of  composites  and  ultimately  be  used  in  op¬ 
timizing  a  composite  drive  shaft  design. 

Section  2  defines  the  scope  and  limitations  of  the 
analysis  and  optimization. 

Section  3  discusses  the  analysis  itself,  the  equations 
and  the  modes  of  failure  that  were  considered.  Examples  of 
the  analysis  are  shown. 

Section  4  briefly  describes  the  optimization  program 
[1,2]  and  presents  the  results  of  optimizing  a  number  of 


shaft  designs.  Two  variable  function  space  diagrams  are 
provided  to  illustrate  the  design  space  for  different  load 
conditions. 

Section  5  contains  recommendations  for  future  investi¬ 
gations  . 

Appendix  A  contains  a  description  of  the  FORTRAN  program 
used  in  the  analysis  and  the  program  text.  Appendix  B  shows 
the  derivation  of  equations  for  deflection,  moment  and  shear 
force . 
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II.  SCOPE  AND  LIMITATIONS 


The  analysis  applies  to  shafts  made  of  isotropic  ma¬ 
terials  with  simply  supported  or  fixed-fixed  end  conditions. 
Basic  shaft  design  formulas  are  used  except  that  the  de¬ 
flection  and  moment  equations  are  derived  to  include  the 
effects  of  a  small  mass  imbalance.  The  shaft  is  assumed 
to  rotate  with  synchronous  whirl. 

To  insure  a  continuous  design  field,  the  optimization 
is  constrained  to  subcritical  speeds  only.  Supercritical 
speeds  introduce  feasible  design  regions  that  are  disjoint 
from  the  primary  design  space.  In  this  investigation,  these 
secondary  areas  are  considered  infeasible. 

Since  the  optimization  tends  to  a  large  radius-small 
thickness  design,  radial  stresses  are  neglected.  However, 
buckling  of  thin  cylinders  due  to  torsion  and  compression 
is  considered.  The  equations  used  for  these  failure  criteria 
do  not  include  dynamic  effects.  A  constant  torque  and  a 
steady,  uniform  axial  force  are  the  loads  considered  in 
the  formulas  [3]. 

The  analysis  incorporates  the  capability  to  design  a 
shaft  that  may  be  used  in  two  or  more  loading  conditions. 

For  example,  a  shaft  may  be  designed  to  transmit  50  HP  at 
300  RPM  with  an  axial  load  of  2000  lbs.  as  well  as  to  transmit 
150  HP  at  3000  RPM  with  a  1000  lb.  axial  load. 
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III.  ANALYSIS 

The  problem  is  a  shaft  required  to  transmit  a  specified 
horsepower  at  a  given  speed.  There  may  be  an  axial  load,  F 
or  an  internal  presence,  P.  (See  Figure  1.)  Basic  equations 
are  used  throughout  the  analysis  and  fundamentals  such  as 
formulas  for  area,  volume,  moment  of  inertia,  polar  moment 
of  inertia  and  radius  of  gyration  are  not  repeated  here. 
However,  certain  stress  calculations  will  be  emphasized  to 
illustrate  the  flow  of  analysis. 

Appendix  B  shows  the  derivations  of  the  equations  for 
bending  moment,  deflection  and  shear  for  pinned-pinned  and 
clamped-clamped  end  conditions.  These  derivations  have 
taken  into  consideration  the  effects  of  a  small  mass  imbalance 
in  conjunction  with  the  rotation  of  a  shaft  with  synchronous 
whirl. 

The  inch-pound-second  system  of  units  (IPS)  is  used  in 
the  computer  program  analysis.  However,  by  changing  the 
equation  for  torque  (line  50  in  the  FORTRAN  text.  Appendix  A) 
to  reflect  the  International  System  of  Units,  the  program 
may  be  used  effectively  with  consistent  SI  units. 
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A.  SIMPLY  SUPPORTED  BEAM 


The  deflected  shape  of  the  beam  is  assumed  as: 

7TX 

6  ~  6  sin  -3-  (see  Fig.  2) .  The  solution  of  the  govern- 
ing  equations  equations  is  given  in  Appendix  B.  The  deflection 
is  maximum  at  x  =  Z/2 


Kxe  +  K2 

[EI-K1(|)4  -  F(i)2] 


(3.1a) 


2 

where:  =  pAgj 

p  =  mass  per  unit  volume 
A  =  cross  sectional  area 
co  =  shaft  speed  is  radians  per  second 
K  *  pAg  =  weight  per  unit  length 

e  =  eccentricity  of  mass  with  respect  to  the  axis 
of  rotation 

Z  =  shaft  length 

E  =  Young's  Modulus 

I  =  Moment  of  Inertia 

F  =  Axial  Load 

Certain  characteristics  of  this  equation  are  of  particular 
interest.  If  there  is  no  rotation  nor  axial  load,  the  formula 
reduces  to  the  classical  equation  for  maximum  beam  deflection. 


(3.1b) 


If  there  is  no  axial  load,  the  equation  becomes  unstable 

2  j?  4 

(deflection  goes  to  infinity)  when  El =  pAw  (£)  or,  solving 
for  the  critical  speed. 


"■ =  TO 


(3.1c) 


which  defines  the  first  fundamental  frequency  of  a  simply 
supported  shaft. 

If  there  is  no  rotation,  instability  arises  when  El  *  F( 


F  * 
c 


2 

tTEI 


(3. Id) 


recognizable  as  Euler’s  column  buckling  criteria. 

With  both  rotation  and  axial  load,  the  critical  speed 


becomes 


u>  =  t r 
c 


2  EI  -  f4>2 


(3 .  le) 


To  avoid  the  instability  regions,  the  denominator  of 
equation  (3.1a)  must  be  given  a  lower  bound  greater  than  zero 
(point  A  in  Fig.  3).  Additionally,  a  maximum  deflection  has 
to  be  imposed  to  prevent  computer  overflow  (point  B  in  Fig.  3) 
These  limits  have  no  relevance  with  the  deflection  constraint 
used  in  the  design  optimization  process  (point  C  in  Fig.  3) , 
but  are  used  only  to  prevent  numerical  ill-conditioning  on 
the  computer. 
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Figure  3 


Limits  on  the  deflection.  The  solid  curve  repre¬ 
sents  the  deflection  values  used  in  the  analysis. 


L6 


The  bending  moment  corresponding  to  the  deflection  calcu¬ 


lated  by  equation  (3.1a)  at  x  -  A/2  is 

K  A^ 

«max  *  K1  1  f  +  S(?>21  *  FS  +  -V-  l3'l£» 

When  the  weight  tends  to  cancel  the  deflection,  the  corre¬ 
sponding  moment  at  the  same  cross  section  becomes 


M  .  »  M 

mm  max 


k2a 


(3 . lg) 


The  shear  force  is  maximum  at  x 


max 


K,  6 A  K.eA  K2A 

— —  +  -=x —  +  1 


( 3 .  lh) 


B.  CLAMPED-CLAMPED  BEAM 

The  deflected  shape  of  the  beam  is  assumed  as: 

<5 

5  »  ~y  [  1  -  cos  — .  (See  Fig.  2b.)  The  deflection  is 
maximum  at  x  *  A/2, 

r  _  «4  r  Kle  +  K2  i  (3.2a) 

°  a  l - - 3 - 5—3 - J 

384  El  -  0.746  K.  A*  -  96  (-)  F 

1  TT 

As  in  the  simply  supported  beam,  similar  observations 
regarding  special  loading  conditions  are  apparent.  If  there 
is  neither  rotation  nor  axial  load,  the  formula  for  maximum 
beam  deflection  results 
k2A 4 

'"max  *  384  El 


(3.2b) 


With  no  axial  load,  instability  occurs  when 
384  El  »  0.746  pAoj2  or 


oo  =  22.68 

c 


I 


El 

pA£4 


(3.2c) 


which  is  the  first  fundamental  frequency  for  a  fixed-fixed 

beam  based  on  the  assumed  deflection  used  here. 

When  there  is  no  rotation,  instability  arises  when 
2 

384  El  =  96  or 

it 


,  _  4tt2EI 

C'  *2  ' 


(3. 2d) 


Euler's  column  buckling  criteria. 

With  rotation  and  axial  load,  the  instability  occurs  at 
a  smaller  speed  than  that  given  by  equation  (3.2c) , 

% 

El  -  0.02533  F&4  (3.2e) 

PA*4 

The  bending  moment  at  x  =  J l/2  corresponding  to  the 
deflection  given  by  equation  (3.2a)  is 

K.  Z  t  n  -j ,  Z 

Mmax  =  IT  I-5-+  —  +e^]+F6+  (3.2f) 

A  lesser  bending  moment  occurs  at  the  same  cross  section 
when  the  weight  tends  to  cancel  the  deflection. 
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M  . 
min 


(3.2g) 


2 

k2jt 

~  ^max  12 

The  maximum  shear  occurs  at  x  =  0. 

K,  Si,  K.ei. 

v  =  -L —  +  -L —  +  -4—  (3. 2h) 

max  4  2  2  v  ' 

The  preceding  calculations  are  performed  in  Subroutine 
BEND,  given  in  Appendix  A. 

C.  STRESS  ANALYSIS 

With  the  deflection,  shear  and  bending  moments  known  at 
the  critical  cross  sections,  the  analysis  proceeds  with  the 
calculation  of  the  stresses  at  the  critical  stress  elements. 
At  midspan,  the  critical  stress  elements  are  on  the  outer 
surface  of  the  shaft  where  the  axial  stress  (tensile  or 
compressive)  is  maximum  (Points  A  and  B,  Fig.  4).  Since  the 
shaft  is  in  synchronous  whirl,  the  same  stress  element,  A, 
will  generally  tend  to  be  in  tension  throughout  the  rotation, 
with  the  magnitude  of  the  axial  stress  varying  as  the  weight 
effect  becomes  alternatingly  additive  and  subtractive. 

Element  A  is  in  the  direction  of  the  mass  eccentricity  from 
the  center  of  the  shaft  while  element  B  is  on  the  outer 
surface  in  the  opposite  direction. 

Most  of  the  following  equations  are  found  in  basic  design 
and  applied  mechanics  books.  References  will  be  noted  when 


necessary. 


Axial  stress  due  to  bending  moment  (tension/compression) 


M  r 

a  (max)  ■  -ma*  — 
x  1 

ox  (min}  -  Mmi;  r° 


(3.3a) 


Axial  stress  due  to  axial  load  (compression) 


a 


x 


L 

A 


(3.3b) 


The  axial  stress  due  to  internal  pressure  for  thin-wall 
cylinders  is: 


o 


X 


Pci 

TT 


(3.3c) 


Superposition  with  careful  regard  to  sign  conventions 
provides  the  maximum  and  minimum  tensile  and  compressive 
axial  stresses. 

Hoop  stress  is  primarily  due  to  internal  pressure  and 
rotational  inertia  [4] 


a 


H 


2  2 
Yw  rQ 


3+\>  [1  +  ( 

4 


l-v 


r .  , 
Hji)2] 
o 


(3. 3e) 


where  Y  *  weight  per  unit  volume 

»  rotation  in  radians  per  second 

v  *  Poisson's  ratio 

r^  »  inside  radius 

r  ■  outside  radius 
o 

g  •  acceleration  due  to  gravity 
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The  shear  stress  at  x  -  1/2  is  caused  only  by  torsion 
(shear  force  is  zero) . 


T 


T 


(3. 3f ) 


where  xT  =  shear  stress  due  to  torsion 
T  =  torque 

J  =  polar  moment  of  inertia 

At  x  =  0,  the  critical  element  is  where  the  shear  stress 
due  to  the  shear  force  along  the  cross  section  is  additive 
with  respect  to  the  shear  stress  due  to  torsion  (element  C, 
Fig.  3) .  The  shear  stress  due  to  Vmax  for  thin  wall  cylinders 
[5]  is 


2  V 


Tv  " 


max 


(3.3g) 


Superposition  of  equations  (3.3f)  and  (3.3g)  gives  the  total 
shear  stress  at  the  critical  stress  element. 

The  axial  stress  at  element  C  due  to  bending  moment  is 
zero.  The  axial  stress  at  this  point  is,  therefore,  caused 
only  by  the  axial  load  and  the  internal  pressure,  equations 
(3.3b)  and  (3.3c). 

The  hoop  stress  is  calculated  by  equations  (3.3d)  and 
(3. 3e)  . 
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stress  element  A 


element  B 


a)  At  midspan 


c)  Typical  stress  element.  Radial 
stresses  are  neglected. 


Figure  4.  Critical  Stress  Elements 


D.  MODES  OF  FAILURE 

Eight  modes  of  failure  are  considered  in  this  investiga¬ 
tion.  The  first  two  evaluate  the  static  strengths  at  the 
critical  stress  elements  using  Von  Mises  criteria  [5].  At 
midspan , 


°mid  £ 


At  x  =  0 


2  2 

a„  a  +o  +  3 
H  x  x  T 


°o  "  j°H  -  °H  0 


x  +  ax  +  3  tt+v 


(3.4a) 


(3.4b) 


Failure  is  assumed  if  the  yield  strength  of  the  material  is 
equal  to  or  less  than  these  Von  Mises  stresses  multiplied  by 
a  safety  factor. 

A  limit  on  the  deflection  is  imposed.  A  common  practice 
is  to  specify  a  maximum  deflection  in  inch  per  foot  of  shaft 
length.  For  example,  as  used  in  this  report. 


5max  “  °-005  inch/foot  length 


Since  the  optimization  may  lead  to  a  thin  cylinder  design, 
buckling  due  to  torsion  and  compression  are  considered.  The 
critical  torque  and  critical  compressive  stress  for  buckling 
of  thin  cylinders  [3]  are; 


T 


c 


3 


TT  VTe 

(l-v"< 


o775 


(3.4c) 
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a 


Et 


(3.4d) 


c  r^3  (1-v2) ' 

Experimental  data  indicate  that  actual  failure  of  the  cylinders 
usually  occurs  below  fifty  percent  of  the  values  calculated 
from  these  equations.  To  accommodate  this  discrepancy  and 
the  fact  that  the  dynamic  effects  are  not  considered  in  the 
derivation  of  these  equations,  comparatively  higher  safety 
factors  must  be  used  for  these  failure  modes. 

The  shaft  is  also  checked  for  column  buckling  using  Johnson's 
or  Euler's  equations  [5]  as  applicable.  The  actual  slenderness 
ratio  is  compared  with  the  slenderness  ratio  where  transition 
occurs  to  determine  which  formulas  apply.  The  calculated 
critical  load,  using  the  appropriate  equations  for  the  end 
condition,  is  compared  to  the  actual  load,  again  incorporating 
a  safety  factor.  Equations  (3. Id)  or  (3. 2d)  are  used  if  the 
slenderness  ratio  falls  in  the  Euler  range.  In  the  Johnson 
region,  the  following  equations  [5]  apply. 


F 

c 

A 


vb  <!>2 


(3 . 4e) 


where  S  =  yield  strength  in  compression 


(-) 

'K' 


slenderness  ratio 
S. 


[— ]2  — 
1 2  rr J  nE 


b 

E  =  Young's  Modulus 
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n  =  1  for  simple  supports 
4  for  fixed-fixed  ends 

The  effects  of  the  fluctuating  axial  stresses  are  evaluated 
for  fatigue  failure  using  the  Goodman  diagram  [5].  Von  Mises 
stresses  for  the  mean  and  alternating  stresses  are  calculated 
as  follows: 


a 

m 


a  +  a  . 
max  mm 

2 


(3.4f) 


a 

a 


a 

max 


or  . 
mm 

1 


(3 . 4g) 


-  "h  +  °m2  +  3  tt2  <3‘4hl 

°a‘  =  °a  (3-4i) 

The  Goodman  diagram  is  drawn  (Fig.  5)  and  the  safety  factor 
is  determined  analytically. 


Figure  5.  Goodman  Diagram 
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r 


SF 


Se 


(3. 4 j ) 


where  Se  =  fully  corrected  endurance  limit 
Su  =  ultimate  strength 

This  safety  factor  is  required  to  be  greater  than  a  specified 
value. 

Lastly,  failure  is  assumed  if  a  specified  percentage  of 
the  critical  speed  as  calculated  by  equations  (3.1e)  or  (3.2e) 
is  less  than  the  shaft  speed. 

E.  ANALYSIS  EXAMPLES 

The  examples  are  chosen  to  illustrate  several  aspects  of 
the  shaft  design  problem.  The  first  is  a  design  which  vio¬ 
lates  almost  all  of  the  constraints  while  the  second  example 
is  one  within  the  feasible  region  of  the  design  field. 

The  material  used  in  all  the  examples  in  this  report  is 
UN  G10150  HR  steel,  the  properties  of  which  are  listed  in 
Table  I  [5] . 
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TABLE  I 


MATERIAL  PROPERTIES 


UNS  G10150  HR  STEEL 


Young ' s  Modulus 
Shear  Modulus 
Poisson's  Ratio 
Weight  Per  Unit  Volume 
Yield  Strength 
Ultimate  Strength 


30. 0E  06  psi 
11. 5E  06  psi 
0.292 

0.282  lbs/in 
27. 0E  03  psi 
50. 0E  03  psi 


1.  Analysis  Example  1 


A  simply  supported  shaft  with  0.1  inch  eccentricity 
is  to  transmit  150  HP  at  3000  RPM  (3150  in-lb  torque) , 
with  an  axial  load  of  2000  lbs.  For  a  shaft  length  of 

120  inches,  inside  radius  of  1.0  inch  and  thickness  of 
1.0  inch,  the  design  violates  the  static  and  dynamic 
strength,  deflection,  buckling  due  to  compression,  and 
the  speed  constraints.  The  analysis  of  this  design  is  given 
below  where  an  asterisk  {*)  denotes  violated  constraints  on 
the  design. 


Volume 

Weight 

Moment  of  Inertia 
Polar  Moment  of  Inertia 
♦Critical  Speed 
Critical  Axial  Load 
Critical  Buckling  Torque 
♦Critical  Buckling  Stress 
♦Deflection 
Bending  Moment 
Shear  Force 
Axial  Stress  at  A 
Axial  Stress  at  B 
Hoop  Stress 

Torsional  Shear  Stress 
Total  Shear  Stress 
Mean  Axial  Stress 
Alternating  Axial  Stress 
♦Von  Mises  at  Midspan 
♦Von  Mises  at  x  =  0 


0.11310E  04  in3 
0.31893E  03  lbs 
0.U781E  02  in* 
0.23562E  02  in4 
0.14767E  04  RPM 
0.18764E  06  lbs 
0.47469E  08  in  lb 
0.905493  07  psi 
0.24000E  02  in 
0.36420E  08  in  lb 
0.62722E  06  lb 
0.95779E  07  psi 
-0.95783E  07  psi 
0.25013E  03  psi 
0.26738E  03  psi 
0.13337E  06  psi 
0.95771E  07  psi 
0.81250E  03  psi 
0.95785E  07  psi 
0.23100E  06  psi 
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The  calculated  deflection  is  the  maximum  value 
imposed  to  avoid  numerical  ill-conditioning  (point  B,  Fig.  3). 
Since  the  deflection  affects  most  of  the  relevant  parameters 
in  the  analysis,  the  calculated  stresses  are  unreliable  and 
not  valid.  These  values,  however,  tend  to  direct  the  design 
towards  the  feasible  region. 

2.  Analysis  Example  2 

This  shaft  is  similar  to  the  first  example  in  design 
conditions  and  loading.  The  only  difference  is  that  the 
dimensions  are  such  that  no  constraints  are  violated,  but 
two  constraints  are  active  or  critical.  The  shaft  is 
120  inches  long  with  an  inside  radius  of  6.39  inches  and 
thickness  of  0.0369  inch.  This,  in  fact,  is  the  optimized 
design  for  the  given  loading  as  calculated  by  COPES-CONMIN 
[1,2].  The  results  of  the'  analysis  follows,  where  a  double 
asterisk  (**)  denotes  values  at  their  design  limit. 


Volume 

0.17866E 

03 

in3 

Weight 

0.50382E 

02 

lbs 

Moment  of  Inertia 

0.30629E 

02 

in4 

Polar  Moment  of  Inertia 

0.61258E 

02 

in4 

♦Critical  Speed 

0.60061E 

04 

RPM 

Critical  Axial  Load 

0.39557E 

05 

lbs 

♦Critical  Buckling  Torque 

0.31490E 

05 

in  lb 

Critical  Buckling  Stress 

0.10400E 

06 

psi 

Deflection 

0. 43813E- 

•01 

in 

Bending  Moment 

0.17026E 

05 

in  lb 

Shear  Force 

0.84895E 

03 

lbs 

Axial  Stress  at  A 

0.22326E 

04 

psi 

Axial  Stress  at  B 

-0.49193E 

04 

psi 

Hoop  Stress 

0.29778E 

04 

psi 

Torsional  Shear  Stress 

0.33079E 

03 

psi 

Total  Shear  Stress 

0.14712E 

04 

psi 

Mean  Axial  Stress 

0.20739E 

04 

psi 

Alternating  Axial  Stress 

0.15872E 

03 

psi 

Von  Mises  at  Midspan 

0.69313E 

04 

psi 

Von  Mises  at  x  =  0 

0.39193E 

04 

psi 
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In  each  of  the  examples,  the  analysis  proceeds  as 
follows.  All  shaft  characteristics  like  cross  sectional 


1 


area,  volume,  weight,  moment  of  inertia  and  polar  moment 
of  inertia  are  first  calculated.  Parameters  that  depend 
only  on  geometry  such  as  critical  torque  and  compressive 
stress  for  buckling  of  thin  cylinders,  and  critical  speed 
are  then  determined  before  Subroutine  BEND,  where  the  de¬ 
flection,  bending  moments  and  shear  force  at  the  critical 
cross  sections  are  computed,  is  called. 

At  midspan,  the  absolute  values  of  the  axial 
stresses  at  elements  A  and  B  (see  Fig.  4a)  are  compared 
to  determine  which  element  is  critical.  In  both  examples, 
element  B  is  considered  critical  for  strength  evaluation. 

At  x  a  0,  the  critical  element,  C,  is  where  the  torsional 
shear  stress  and  the  shear  stress  due  to  the  shear  force 
along  the  cross  section  are  additive.  The  Von  Mises 
stresses  for  each  of  the  critical  elements  are  then  calcu¬ 
lated  using  equations  (3.4a)  and  (3.4b).  The  rest  of  the 
analysis  is  a  straight  forward  application  of  the  equations 
for  the  critical  parameters  and  the  constraints. 
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IV.  OPTIMIZATION,  RESULTS  AND  CONCLUSIONS 

A.  OPTIMIZATION  AND  RESULTS 

There  are  a  number  of  optimization  programs  available, 
each  one  using  different  techniques  in  locating  the  desired 
optimum  design.  COPES-CONMIN  [1,2]  is  a  versatile  program 
that  may  be  used  for  sensitivity  analysis  and  two  variable 
function  space  study  as  well  as  an  optimization  tool. 

Provided  with  a  user  supplied  analysis  program  (Subroutine 
ANALIZ)  where  the  objective  function,  constraints  and  other 
relevant  parameters  are  calculated,  it  determines  a  usable 
and  feasible  sector  from  which  a  search  direction  is  chosen. 
This  choice  of  search  direction  is  made  by  using  such  informa¬ 
tion  as  the  gradients  of  the  objective  function,  gradients  of 
active  and  violated  constraints  and  the  degree  of  push-off 
these  constraints  violations  generate.  This  iterative  process 
of  minimizing/maximizing  the  objective  function  by  changing 
the  design  variables  is  terminated  when  no  further  improvement 
can  be  made. 

In  this  exercise,  the  objective  function  to  be  minimized 
is  the  volume,  with  the  radius  and  thickness  as  design 
variables.  Several  design  examples  are  shown  to  demonstrate 
the  effects  of  the  different  parameters  on  the  optimized 
design  and  on  the  design  space  itself.  Seven  examples  are 


31 


presented  corresponding  to  the  loading  conditions  as  listed 
in  Table  II.  Two  variable  function  space  diagrams  are 
generated  to  illustrate  these  effects.  The  numbered  curves 
in  the  graphs  represent  the  corresponding  constraints  shown 
in  Table  III.  The  shaded  side  of  a  curve  represents  the  un¬ 
feasible  area  of  the  design  space  with  respect  to  the  cor¬ 
responding  constraint.  The  dashed  curves  represent  contours 
of  the  objective  function,  the  numbers  indicating  the 
volume  in  cubic  inches. 


TABLE  II 

OPTIMIZATION  EXAMPLES 


Loads 


NO. 

Horsepower 

RPM 

Pressure 

(psi) 

Axial  Load 
(lbs) 

1 

150 

300 

0 

0 

2 

300 

1000 

0 

3 

150 

300 

0 

2000 

4 

150 

3000 

0 

0 

5 

150 

3000 

1000 

0 

6 

150 

3000 

0 

2000 

7* 

150 

300 

0 

0 

150 

3000 

0 

2000 

♦Example  No.  7  is  a  design  for  two  load  conditions. 


OO^J<TvUl»UtorOM 


TABLE  III 


Number 


DESIGN  CONSTRAINTS 

Constraint 


Static  Strength  at  Midspan 
Static  Strength  at  x  =  0 
Deflection 

Buckling  due  to  Torque 
Buckling  due  to  Compression 
Column  Buckling 
Dynamic  Strength  (Fatigue) 
Speed 
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EXAMPLE  1.  i  tNlTIAL  DESIGN 
SHAFT  DI MENSI ONS  s 

RI  «  0.10003E  01  TH  »  0.10000E  01  SL  -  0.12000E 
DESIGN  CONDITIONS  : 

ECCENTRICITY  »  0.10000=  30 
END  CONDITIONS  :  P INNED-PT NNED 
LOAOS  i 

HORSEPOWER  RPM  PRESSURE  AXIAL  L3A0 

0.15000E  03  0.30000E  03  0.0  0.0 


ANALYSIS  RESULTS 

T0R3Ue(l) 

VOLUME 

weight 

MOMENT  OF  INERTIA 
POLAR  MOMENT  OF  INERTIA 
CRITICAL  SPEED 
CRITICAL  AXIAL  LOAD 
CRITICAL  BUCKLING  TORQJS 
CRITICAL  BUCCING  STRESS 
DEFLECTION 
SENDING  MOMENT 
SHEAR  FORCE 
AXIAL  STRESS  AT  A 
AXIAL  STRESS  AT  B 
HOOP  STRESS 
TORSIONAL  SHEAR  STRESS 
TOTAL  SHEAR  STRESS 
MEAN  AXIAL  STRESS 
ALTERNATING  AXIAL  STRESS 
VON  MISES  AT  MIDSPAN 
VON  MISES  AT  X*0. 


0.31500E  05 
0.11310E  04 
0.31893F  03 
0.1179  IE  02 
0.23562E  02 
0. 1482 BE  04 
0.18754E  06 
0.47469E  08 
0. 90549E  07 
0.26584E-01 
0.50477E  04 
0.2371 4E  03 
0.85692E  03 
-0.85692S  03 
0.2501 3 E  01 
0.26738E  04 
0.27178E  04 
0.44763E  02 
0.8L216E  03 
0.47095E  04 
0.47073E  04 
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EXAMPLE  i.  :  OPTIMUM  DESIGN 
SHAFT  DIMENSIONS  : 

RI  »  0. 25847E  01  TH  »  0.U124E  00  SL  =  0.12000E 
DESIGN  CONDITIONS  : 

ECCENTRICITY  *  0.10000E  00 
END  CONDITIONS  s  P I NNE3-PI NNEO 
LOADS  : 


HORSEPOWER  RPM 

PRESSURE 

AXIAL  LOAD 

0.150005  03  0.30000E  03 

0.0 

o 

. 

o 

ANALYSIS  RESULTS 

TORQUE <11 

0.31500E 

05 

VOLUME 

0.2214S6 

03 

WEIGHT 

0.62451E 

02 

MOMENT  OF  INERTIA 

0.64353E 

01 

POLAR  MOMENT  OF  INERTIA 

0. 123  71 E 

02 

CRITICAL  SPEED 

0.247576 

04 

CRITICAL  AXIAL  LOAD 

0.45137E 

05 

CRITICAL  BULLING  TORQUE 

0.315006 

06 

CRITICAL  BUCKLING  STRESS 

0.74729E 

06 

DEFLECTION 

0. 92753E- 

-02 

BENDING  MOMENT 

0.95479E 

03 

SHEAR  FORCE 

0. 3963  IE 

02 

AXIAL  STRESS  AT  A 

0.39998E 

03 

AXIAL  STRESS  AT  B 

-0.399996 

03 

HOOP  STRESS 

0.51655E 

01 

TORSIONAL  SHEAR  STRESS 

0.659806 

04 

TOTAL  SHEAR  STRESS 

0.6641 OE 

04 

MEAN  AXIAL  STRESS 

0.75507E 

01 

alternating  axial  stress 

0.39243E 

03 

VON  MISES  AT  MIDSPAN 

0.114356 

05 

VON  MISES  AT  X«0. 

0.11503E 

05 
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EXAMPLE  2.  :  INITIAL  DESIGN 
SHAFT  DIMENSIONS  : 

RI  =  0.10000=  01  TH  =  0.10000E  01  $L  *  0.12000E 
DESIGN  CONDITIONS  : 

ECCENTRICITY  =  0.10000E  00 
END  CONDITIONS  :  PI NNED-PI NNED 
LOADS  : 


HORSEPOWER  RPM 

PRESSURE 

AXIAL  LOAD 

0.15000E  03  0.30000E  03 

0. 10000E  04 

o 

. 

o 

ANALYSIS  RESULTS 

TORQJE(l) 

0.31500E 

05 

VOLUME 

0.  1131 OE 

04 

WEIGHT 

0.31893E 

03 

MOMENT  OF  INERTIA 

0. 1178  IE 

02 

POLAR  MOMENT  OF  INERTIA 

0 . 23562  E 

02 

CRITICAL  SPEED 

0.14828E 

04 

CRITICAL  AXIAL  LOAD 

0.18764E 

06 

CRITICAL  3UC<L I NG  TORQUE 

0.47469E 

08 

CRITICAL  BUCRLING  STRESS 

0.90549E 

07 

DEFLECTION 

0.26584E-01 

BENDING  MOMENT 

0.50477E 

04 

SHEAR  FORCE 

0.2071 4E 

03 

AXIAL  STRESS  AT  A 

0.18569E 

04 

AXIAL  STRESS  AT  B 

0.14309E 

03 

HOOP  STRESS 

0.66917E 

03 

TORSIONAL  SHEAR  STRESS 

0.267385 

04 

TOTAL  SHEAR  STRESS 

0.27178E 

04 

MEAN  AXIAL  STRESS 

0.10448E 

04 

ALTERNATING  AXIAL  STRESS 

0.81216E 

03 

VON  MISES  AT  MIDSPAN 

0.49093E 

04 

VON  MISES  AT  X»0. 

0.47546E 

04 
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0°  TI MUM  DESIGN 


EX\MPLE  2.  s 

SHAFT  DIMENSIONS  : 

RI  *  0.18625E  01  TH  =  0.2S560E  00  SL  =  0.12000E 

DESIGN  CONDITIONS  i 
ECCENTRICITY  *  0.10000E  00 
END  CONDITIONS  :  Pi NNED-PI NNED 
LOADS  J 

HORSEPOWER  RPM  PRESSURE  AXIAL  LOAD 

0.15000E  03  0.30000E  03  0.10000E  04  0.0 


ANALYSIS  RESULTS 


TORQUE ( I ) 

0.31500E 

05 

VOLUME 

0.39958E 

03 

WEIGHT 

0.I1268E 

03 

MOMENT  OF  INERTIA 

0.6S573E 

01 

POLAR  MOMENT  OF  INERTIA 

0.I33ISE 

02 

CRITICAL  S  PEE  0 

0.18754E 

04 

CRITICAL  AXIAL  LOAD 

0.75144E 

05 

CRITICAL  BUCXLING  TOROJE 

0.23553E 

07 

CRITICAL  BUCXLING  STRESS 

0.22602E 

07 

DEFLECTION 

0.I6358E-0I 

BENDING  MOMENT 

0.L7476E 

04 

SHEAR  FORCE 

0.72247E 

02 

AXIAL  STRESS  AT  A 

0.75709E 

04 

AXIAL  STRESS  AT  8 

0.64537E 

04 

HOOP  STRESS 

0.65487E 

04 

TORSIONAL  SHEAR  STRESS 

0.50344E 

04 

TOTAL  SHEAR  STRESS 

0. 50779  E 

04 

MEAN  AXIAL  STRESS 

0.70306E 

04 

ALTERNATING  AXIAL  STRESS 

0.54027E 

03 

VON  MISES  AT  MIOSPAN 

0.U254E 

05 

VON  MISES  AT  X*0 • 

0.109&5E 

05 
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EXAMPLE  3.  :  MITUL  DESIGN 
SHAFT  DIMENSIONS  : 

RI  =  0.100005  01  TH  =  0.10000E  01  SL  *  0.12000E 
DESIGN  CONDITIONS  : 

ECCENTRICITY  =  0.10000E  00 
END  CONDITIONS  :  PI NNED-PI NNED 
LOADS  : 

HORSEPOWER  RPM  PRESSURE  AXIAL  LOAD 

0.15000E  03  0.30000E  03  0.0  0.20000E  04 


ANALYSIS  RESULTS 


TORQUE (11 

0.31500E 

05 

VOLUME 

0.11310E 

04 

WEIGHT 

0.31893E 

03 

MOMENT  OF  INERTIA 

0.1179LE 

02 

POLAR  MOMENT  OF  INERTIA 

0 .235625 

02 

CRITICAL  SPEED 

0.14767E 

04 

CRITICAL  AXIAL  LOAO 

0.13764E 

06 

CRITICAL  BUC< L I NG  TORQUE 

0.474695 

08 

CRITICAL  BUC<1 ING  STRESS 

0.90549E 

07 

DEFLECTION 

0.26814E-01 

BENDING  MOMEU 

0.541445 

04 

SHEAR  FORCE 

0.20720E 

03 

AXIAL  STRESS  AT  A 

0. 70697  E 

03 

AXIAL  STRESS  \T  8 

-0.11314E 

04 

HOOP  STRESS 

0. 2501 3  E 

01 

TORSIONAL  SHEAR  STRESS 

0.26739E 

04 

TOTAL  SHEAR  STRESS 

0.27178E 

04 

MEAN  AXIAL  STRESS 

-0.10519E 

03 

ALTERNATING  AXIAL  STRESS 

0.81216E 

03 

VON  MISES  AT  MIOSPAN 

0. 47577 E 

04 

VON  MISES  AT  X=0 . 

0.47073E 

04 
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EXAMPLE  3 


OPTIMUM  DESIGN 


SHAFT  DIMENSIONS  : 

RI  *  0.25907E  01  TH  *  0. 11119E  00  SL  =  0. 12000E 
OESIGN  CONDITIONS  : 

ECCENTRICITY  *  0.10000E  00 
END  CONDITIONS  :  PI NNED-PI NNED 
LOADS  : 

HORSEPOWER  RPM  PRESSURE  AXIAL  LOAD 

O.ISOOOE  03  0.30000E  03  0.0  0.20000E  04 


ANALYSIS  RESULTS 


TORQUE (1 I 

0.31500E 

05 

VOLUME 

0.22I86E 

03 

WEIGHT 

0.6256VE 

02 

MOMENT  OF  INERTIA 

0.64766E 

01 

POLAR  MOMENT  OF  INERTIA 

0.12953E 

02 

CRITICAL  SPEED 

0.2463SE 

04 

CRITICAL  AXIAL  LOAD 

0.45241E 

05 

CRITICAL  BUC<L I NG  TORQJE 

0.315005 

06 

CRITICAL  BUCKLING  STRESS 

0.74527E 

06 

deflection 

0.93752E-02 

BENDING  MOMENT 

0.9817LE 

03 

SHEAR  FORCE 

0.39758E 

02 

AXIAL  STRESS  AT  A 

-0. 6722 1 E 

03 

AXIAL  STRESS  AT  B 

-0. 1491 3  6 

04 

HOOP  STRESS 

0.51898E 

01 

TORSIONAL  SHEAR  STRESS 

0.6570&E 

04 

TOTAL  SHEAR  STRESS 

0.66136E 

04 

MEAN  AXIAL  STRESS 

-0.10637E 

04 

ALTERNATING  AXIAL  STRESS 

0.3915LE 

03 

VON  MISES  AT  MIDSPAN 

0.11478E 

05 

VON  MISES  AT  X*0. 

0.I1455E 

05 
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INSIDE  RRDIUS  INCHES  ( 


INITIAL  DESIGN 


EXAMPLE  4.  : 

SHAFT  DIMENSIONS  J 
RI  a  0.10000=  01  TH  *  0.10000E  01 
DESIGN  CONDITIONS  : 

ECCENTRICITY  *  0.10000E  00 
END  CONDITIONS  :  PI NNED-PI NNED 
LOADS  : 

HORSEPOWER  RPW  PRESSURE 

0.15000E  03  0.30000E  04  0.0 


SL  «  0.12000E 


AXIAL  LOAD 

0.0 


TORQUE ( 1 ) 

0.31500E 

04 

VOLUME 

0.1131 OE 

04 

WEIGHT 

0.31893E 

03 

MOMENT  OF  INERTIA 

0.11781E 

02 

POLAR  MOMENT  OF  INERTIA 

0.23562E 

02 

CRITICAL  SPEED 

0.14828E 

04 

CRITICAL  AXIAL  LOAD 

0.18764E 

06 

CRITICAL  BUC< u I NG  TORQUE 

0.474S9E 

08 

CRITICAL  BUC<L!NG  STRESS 

0.90549E 

07 

deflection 

0.24000E 

02 

BENDING  MOMENT 

0.23801E 

08 

SHEAR  FORCE 

0.62722E 

06 

AXIAL  STRESS  AT  A 

0.40406E 

07 

AXIAL  STRESS  AT  B 

-0, 40 406 E 

07 

HOOP  STRESS 

0.25013E 

03 

TORSIONAL  SHEAR  STRESS 

0.2&738E 

03 

total  SHEAR  STRESS 

0.13337E 

06 

MEAN  AXIAL  STRESS 

0.40397E 

07 

ALTERNATING  AXIAL  STRESS 

0.81200E 

03 

VON  MISES  AT  MIDSPAN 

0.40404E 

07 

VON  MISES  AT  X*0. 

0.23100E 

06 

i 

y 

i 

t 

i 
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EXAMPLE  4 


OPTIMUM  DESIGN 


SHAFT  DIMENSIONS  : 

RI  *  0.63798E  01  TH  *  0.36963E-01  SL  =  0.12000E 
DESIGN  CONDITIONS  : 

ECCENTRICITY  »  0.100005  00 
END  CONDITIONS  :  PINNED-PI NNED 
LOADS  : 

HORSEPOWER  RPM  PRESSURE  AXIAL  LOAD 

0.15000E  03  0.30000E  04  0.0  0.0 


ANALYSIS  RESULTS 


TORQUE (1) 

0.31500E 

04 

VOLUME 

0.17831E 

03 

WEIGHT 

0.50283E 

02 

MOMENT  OF  INERTIA 

0.3041 5E 

02 

POLAR  MOMENT  OF  INERTIA 

0 .6083  3  E 

02 

critical  SPEED 

0.60005  E 

04 

CRITICAL  AXIA.  LOAD 

0.39476E 

05 

CRITICAL  8UCXLING  TORQUE 

0.31494E 

05 

critical  buc<ling  STRESS 

0.10432E 

06 

DEFLECTION 

0.4392  2  E-01 

BENDING  MOMENT 

0.76214E 

04 

SHEAR  FORCE 

0.84773E 

03 

AXIAL  STRESS  AT  A 

0. 1607 9 E 

04 

AXIAL  STRESS  AT  B 

-0.16079E 

04 

HOOP  STRESS 

0.29628E 

04 

TORSIONAL  SHEAR  STRESS 

0.33228E 

03 

TOTAL  SHEAR  STRESS 

0.1473  3E 

04 

MEAN  AXIAL  STRESS 

0.14488E 

04 

ALTERNATING  AXIAL  STRESS 

0.1591 3E 

03 

VON  MISES  AT  MIDSPAN 

0.26327E 

04 

VON  MISES  AT  X*0. 

0.39103E 

04 
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INSIDE  RADIUS  INCHES 


EXAMPLE  5.  : 

INITIAL  DESIGN 

SHAFT  DIMENSIONS  : 

RI  *  0. 100005  01  TH  * 

O.IOOOOE  01  SL 

*  0.12000E 

DESIGN  CONDITIONS  : 

ECCENTRICITY  *  O.IOOOOE 

00 

END  CONDITIONS  :  PINNED- 

PINNED 

LOADS  J 

HORSEPOWER  RPM 

PRESSURE 

AXIAL  LOAD 

0.15000E  03  0.30000E  04  O.IOOOOE  04 

0.0 

ANALYSIS  RESULTS 

TOROUE ( 1 ) 

0.31500E 

04 

VOLUME 

0. 1131 OE 

04 

WEIGHT 

0.31893E 

03 

MOMENT  OF  INERTIA 

0.11731E 

02 

POLAR  MOMENT  OF  INERTIA 

0.23562E 

02 

CRITICAL  SPEED 

0.14828E 

04 

CRITICAL  AXIAL  LOAD 

0.18764E 

06 

CRITICAL  8UC<L ING  TORQUE 

0.47469E 

08 

CRITICAL  BUCKLING  STRESS 

0.90549E 

07 

DEFLECTION 

0.24000E 

02 

BENDING  MOMENT 

0.23801 E 

08 

SHEAR  FORCE 

0.6272  2 E 

06 

AXIAL  STRESS  AT  A 

0.4041 6E 

07 

AXIAL  STRESS  AT  B 

-0. 40396E 

07 

HOOP  STRESS 

0.91680E 

03 

TORSIONAL  SHEAR  STRESS 

0.26738E 

03 

TOTAL  SHEAR  STRESS 

0.13337E 

06 

MEAN  AXIAL  STRESS 

0.40407E 

07 

ALTERNATING  AXIAL  STRESS 

0.81200E 

03 

VON  USES  AT  MIDSPAN 

0.4041 l E 

07 

VON  MISES  AT  X*0. 

0.23100E 

06 
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EXAMPLE  5.  :  OPTIMUM  DESIGN 
SHAFT  DIMENSIONS  : 

RI  «  0.60844E  01  TH  *  0.61266E  00  SL  =  0. 12000E 
DESIGN  CONDITIONS  : 

ECCENTRICITY  *  0.10000E  00 
END  CONDITIONS  :  P INNEO-PI NNEO 
LOADS  : 


HORSEPOWER  RPM 

PRESSURE 

AX 

0.15D00E  03  0.30000E  04 

0.10000E  04 

O.i 

ANALYSIS  RES UL 

TS 

TORQUE (1 ) 

0.31500E 

04 

VOLUME 

0. 2952  IE 

04 

WEIGHT 

0. 8325 OF 

03 

MOMENT  OF  INERTIA 

0.50353E 

03 

POLAR  MOMENT  OF  INERTIA 

0.10071E 

04 

CRITICAL  SPEED 

0.60003E 

04 

CRITICAL  AXIAL  LOAD 

0.65357E 

06 

CRITICAL  BUCKLING  T3RQJE 

0.34401E 

08 

CRITICAL  BUCKLING  STRESS 

0.16557E 

07 

DEFLECTION 

0.43926E- 

■01 

BENDING  MOMENT 

0.1261 9E 

06 

SHEAR  FORCE 

0.14035E 

05 

AXIAL  STRESS  AT  A 

0.11609E 

05 

AXIAL  STRESS  AT  B 

0.82526E 

04 

HOOP  STRESS 

0.12589E 

05 

TORSIONAL  SHEAR  STRESS 

0.20948E 

02 

TOTAL  SHEAR  STRESS 

0. 1162  OE 

04 

MEAN  AXIAL  STRESS 

0.11443E 

05 

ALTERNATING  AXIAL  STRESS 

0.16609E 

03 

VON  MISES  AT  MIDSPAN 

0.12129E 

05 

VON  MI SES  AT  X*0. 

0.12749E 

05 

i 

i 


INSIDE  RADIUS 


EXAMPLE  6.  :  INITIAL  DESIGN 
SHAFT  DIMENSIONS  i 

R I  *  0.13000:  01  TH  *  0. 10000E  01  SL  ■  0.12000E 
OESIGN  CONDITIONS  : 

ECCENTRICITY  «  0.10000E  00 
END  CONDITIONS  :  PI NNED-PI NNED 
LOADS  : 


HORSEPOWER  RPM 

PRESSURE 

AXIAL  LOAD 

0.15000E  03  0.30000E  04 

0.0 

0.20003E  04 

ANALYSI S  RESULTS 

TORQJE ( 1 J 

0.31500E 

04 

VOLUME 

0. 1131 OE 

04 

WEIGHT 

0.31893E 

03 

MOMENT  OF  INERTIA 

0. 11 73  IE 

02 

POLAR  MOMENT  OF  INERTIA 

0.23552E 

02 

CRITICAL  SPEED 

0.14757E 

04 

CRITICAL  AXIA.  LOAO 

0.18764E 

06 

CRITICAL  BUC<LING  TORQUE 

0.47469E 

08 

CRITICAL  8UCXLING  STRESS 

0.90549E 

07 

deflection 

0.24000E 

02 

BENDING  MOMENT 

0. 5642  OE 

08 

SHEAR  FORCE 

0.62722E 

06 

AXIAL  STRESS  AT  A 

0.95779E 

07 

AXIAL  STRESS  AT  B 

-0.95783E 

07 

HOOP  STRESS 

0.25013E 

03 

TORSIONAL  SHEAR  STRESS 

0.26733E 

03 

TOTAL  SHEAR  STRESS 

0.13337E 

06 

MEAN  AXIAL  STRESS 

0.95771E 

07 

ALTERNATING  AXIAL  STRESS 

0.81250E 

03 

VON  MI SES  AT  MIDSPAN 

0. 95785  E 

07 

VON  MI SES  AT  X*0. 

0.231005 

06 
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EXAMPLE  6 


OPTIMUM  DESIGN 


SHAFT  DIMENSIONS  : 

RI  »  0.63959E  01  TH  a  D.36942E-01  SL  *  0.12000E 
DESIGN  CONDITIONS  : 

ECCENTRICITY  «  0.100005  00 
END  CONDITIONS  :  PI NNED-PI NNEO 
LOADS  s 


HORSEPOWER  RPM 

PRESSURE 

AXIAL  LOAD 

0.15000E  03  0.30000E  04. 

0.0 

0.20000E  04 

ANALYSIS  RESULTS 

TOROJE (1  ) 

0.31500E 

04 

VOLJME 

0.1*866= 

03 

WEIGHT 

0.503825 

02 

MOMENT  OF  IN  FR  T I  A 

0.306295 

02 

POLAR  MOMENT  OF  INERTIA 

0 .6 1253  E 

02 

CRITICAL  SPEED 

0.6006  IE 

04 

CRITICAL  AXIAL  LOAD 

0.39557E 

05 

CRITICAL  BUCXLTNG  TORQUE 

0.31490E 

05 

CRITICAL  BUCXLING  STRESS 

0.10400E 

06 

DEFLECTION 

0.43813E- 

-01 

BENDING  MOMENT 

0.17026E 

05 

SHEAR  FORCE 

0.84895= 

03 

AXIAL  STRESS  AT  A 

0.22325E 

04 

AXIAL  STRESS  AT  B 

-0.49193E 

04 

HOOP  STRESS 

0.29773E 

04 

TORSIONAL  SHEAR  STRESS 

0.33079E 

03 

TOTAL  SHEAR  STRESS 

0.14712E 

04 

MEAN  AXIAL  STRESS 

0.20739E 

04 

ALTERNATING  AXIAL  STRESS 

0.15372E 

03 

VON  MISES  AT  MIDSPAN 

0.69313E 

04 

VON  MISES  AT  X*0. 

0.39193E 

04 
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EXAMPLE  7 


INITIAL  DESIST 


SHAFT  DIMENSI DNS  t 

RI  »  0.10000=  01  TH  =  0.10000E  01  SL  =  0.12000E 
DESIGN  CONDITIONS  : 

eccentricity  =  0.10000=  00 

END  CONDITIONS  :  P I NNED-P I NN ED 
LOADS  s 

HORSEPOWER  RPM  PRESSURE  AXIAL  LOAD 

0.15000E  03  0.30000E  03  0.0  0.0 

0.15000E  03  0.30000E  04  0.0  0.20000E  04 


ANALYSIS  RESULTS 


TORQUE  UJ 

0.31500E 

05 

TORQUE (2 ) 

0.31500E 

04 

VOLUME 

0 . 1131 0  E 

04 

WEIGHT 

0.31893E 

03 

MOMENT  OF  INERTIA 

0.117B1E 

02 

POLAR  MOMENT  OF  INERTI A 

0.23562E 

02 

CRITICAL  SPEED 

0.14767E 

04 

CRITICAL  AXIAL  LOAD 

0. 137S4E 

06 

CRITICAL  BUCKLING  TORQUE 

0.47469E 

08 

CRITICAL  BUCXL I NG  STRESS 

0.90549E 

07 

DEFLECTION 

0.24000E 

02 

BENDING  MOMENT 

0.56420E 

08 

SHEAR  FORCE 

0.62722= 

06 

AXIAL  STRESS  AT  A 

0.95779E 

07 

AXIAL  STRESS  AT  B 

-0.95783E 

07 

HOOP  STRESS 

0.25013E 

03 

TORSIONAL  SHEAR  STRESS 

0.26733E 

03 

TOTAL  SHEAR  STRESS 

0.13337= 

06 

MEAN  AXIAL  STRESS 

0.95771= 

07 

ALTERNATING  AXIAL  STRESS 

0.81250E 

03 

VON  MI  SES  AT  MIOSPAN 

0.95785E 

07 

VON  MI SES  AT  X*0. 

0.23100E 

06 
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EXAMPLE  7.  j  DPTIMUM  DESIGN 
SHAFT  DIMENSIONS  i 

RI  =  0. 63 55  IE  01  TH  =  0.93009E-01  SL  =  0.12000E 
DESIGN  CONDITIONS  : 

ECCENTRICITY  *  0.10000=  30 
END  CONDITIONS  :  P INNED-PI NNED 
LOADS  : 


HORSEPOWER  RPM 

PRESSURE 

AXIAL  LOAO 

0.15000E  03  0.30000E  03 

0.0 

o 

• 

o 

0.150005  03  3.300006  04 

0.0 

0. 20D0DE  04 

ANALYSIS  RESULTS 

TORQUE ( 1 )  0.31533E 

05 

TORQUE (2) 

0 .315005 

04 

VOLUME 

0.44893E 

03 

WEIGHT 

0 •  12663E 

03 

MOMENT  OF  INERTIA 

0.76653E 

02 

°0L AR  MOMENT  OF  INERTIA 

0.153325 

03 

CRITICAL  SPEED 

0.60003E 

04 

CRITICAL  AXIA.  LOAD 

0.99393E 

05 

CRITICAL  BUCXLING  TORQUE 

0. 31 571 E 

06 

CRITICAL  BUCXLING  STRESS 

0.261226 

06 

deflection 

0.43932E- 

•01 

BENDING  MOMENT 

0.428935 

05 

SHEAR  FORCE 

0.21344E 

04 

AXIAL  STRESS  AT  A 

0.3073  3  E 

04 

AXIAL  STRESS  AT  B 

-0.41425E 

04 

HOOP  STRESS 

0.29828E 

04 

TORSIONAL  SHEAR  STRESS 

0.13243E 

03 

TOTAL  SHEAR  STRESS 

0.12736E 

04 

MEAN  AXIAL  STRESS 

0.29136E 

04 

alternating  axial  stress 

0.15973E 

03 

VON  MI SES  AT  MIDSPAN 

0.620215 

04 

VON  MI SES  AT  X*0. 

0.37099E 

04 
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B.  OBSERVATIONS  AND  CONCLUSIONS 


With  all  other  parameters  held  constant,  an  increase  in 
speed  drives  the  design  to  a  very  large  inside  radius  while 
reducing  the  thickness  considerably.  This  trend  is  so 
strong  that  a  reasonable  lower  bound  on  the  thickness  must 
be  given  to  prevent  numerical  ill-conditioning.  The  intro¬ 
duction  of  eccentricity  increases  the  effect  of  rotation  on 
the  stresses  and  the  optimization  produces  slightly  larger 
values  for  the  design  variables.  This  effect  is  reduced, 
however,  by  the  constraint  on  the  deflection. 

As  the  internal  pressure  is  increased,  the  tendency  is 
to  produce  an  optimized  radius  that  is  smaller,  and  a  larger 
thickness.  This  may  be  explained  by  the  direct  effect  of 
the  internal  pressure  on  the  hoop  stress  which  in  turn  would 
require  a  larger  thickness  to  support  the  higher  stress. 

The  effect  of  axial  load  is  reflected  on  the  deflection. 
Again,  if  very  small  deflections  are  allowed,  the  axial  load 
effects  are  reduced.  The  major  parameter  that  is  affected 
is  the  axial  stress. 

A  notable  characteristic  shown  by  the  two  variable 
function  space  diagrams  is  the  tendency  of  the  strength 
constraints  to  curl  back  as  the  radius  is  increased.  This 
signifies  the  existence  of  an  upper  limit  on  size  for  a  given 
design  load.  The  increase  in  mass,  especially  at  high  speeds, 
increases  the  stresses  considerably  leading  to  violation  of 
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the  strength  constraints.  This  phenomenon  is  primarily  a 
result  of  the  eccentricity  of  the  mass  but  also  is  due  to 
the  static  weight  of  the  shaft  itself. 

It  is  also  noted  that  these  strength  constraints  generally 
run  close  to  each  other  in  an  almost  parallel  manner.  The 
dynamic  or  fatigue  constraint  maintains  a  position  such  that 
coming  from  the  feasible  range  towards  the  infeasible  area, 
it  is  activated  before  the  static  constraints  are  encountered. 
In  effect,  the  latter  constraints  are  redundant  in  rotating 
shafts.  However,  these  are  kept  in  the  analysis  to  allow 
solution  of  static  problems. 

The  versatility  of  the  program  should  not  be  lost  in  the 
emphasis  on  drive  shaft  applications.  With  the  appropriate 
loadings,  one  may  use  it  to  solve  structural  problems  like 
columns,  beams  and  pressure  vessel  designs.  The  objective 
function  need  not  be  confined  to  the  volume.  For  example, 
given  a  specific  initial  dimensions,  one  may  minimize  the 
Young's  modulus  to  determine  whether  or  not  a  certain 
material  may  be  utilized;  or  to  determine  how  low  the  modulus 
could  get  before  the  shaft  fails.  If  the  outer  dimension 
is  critical  to  the  design,  one  may  set  a  specific  value  on 
the  diameter  and  let  all  other  parameters  vary  during 
optimization. 

In  summary,  numerical  optimization  provides  an  efficient 
and  effective  way  of  improving  shaft  designs.  The  analysis 
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used  in  this  exercise  is  a  first  step  when  compared  to  the 
boundless  refinements  in  techniques  that  are  possible. 


V.  RECOMMENDATIONS  for  future  investigations 

The  study  has  shown  the  feasibility  of  using  numerical 
optimization  techniques  in  the  design  of  drive  shafts  within 
the  limitations  imposed  in  the  analysis.  Further  studies  on 
the  same  design  field  may  be  pursued  by  eliminating  some  of 
these  limitations.  For  example,  the  unfeasible  region  defined 
by  the  speed  constraint  extends  only  to  some  discrete  distance 
and  a  supercritical  feasible  design  area  does  exist.  It  will 
be  interesting  to  develop  optimization  techniques  to  accommo¬ 
date  this  disjoint  design  field  and  to  investigate  the  pecu¬ 
liarities  this  will  introduce  to  the  analysis. 

Another  investigation  may  be  concerned  with  shafts  with 
variable  thickness  and  radius  along  its  length.  The  inclusion 
of  radial  stresses  in  the  calculations,  the  use  of  finite 
element  methods  in  three  dimensional  stress  analysis  of  the 
problem,  and  vibration  considerations  are  but  some  of  the 
aspects  of  research  that  may  be  useful  in  future  analyses  of 
composite  materials. 

The  challenge  of  including  the  unique  material  character¬ 
istics  of  composites  in  the  design  will  invariably  stimulate 
a  lot  of  young  minds.  However,  the  handling  of  these 
direction  oriented  properties  in  stress  analysis  have  already 
been  investigated  [6  through  12],  and  the  tools  required  to 
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carry  out  the  program  are  available.  COMAND  [13]  offers  a 
simplified  composite  analysis  procedure  for  panels.  Such 
procedures,  adapted  with  the  program  used  in  this  investi¬ 
gation  could  result  in  a  simplified  approach  to  composite 
drive  shaft  analysis.  In  this  procedure,  failure  of  the 
shaft  is  assumed  if  a  single  ply  fails  and  the  entire 
analysis  is  performed  in  a  ply  by  ply  basis.  Since  the 
stress  elements  have  to  be  in  a  differential  level,  the 
composite  material  is  assumed  to  be  homogenous  and  uni- 
directionally  isotropic  at  the  differential  range.  It  will 
also  be  necessary  to  assume  that  stresses  throughout  the 
thickness  of  the  ply  is  uniform  and  equal  to  the  stresses 
at  the  outer  element  of  the  ply.  Analysis  may  then  be 
performed  in  the  classical  procedure,  solving  for  the 
stresses  along  the  major  shaft  axes.  The  material  properties 
to  be  used  in  the  calculation  must  be  equivalent  values 
derived  from  the  directional  properties  using  transformation 
of  axes  equations.  References  10  and  13  present  procedures 
to  calculate  the  corresponding  material  properties,  stresses 
and  strains  along  the  major  axes  of  the  shaft  and  the  fiber 
orientation  coordinates. 

More  involved  investigations  could  deal  with  failure 
modes  of  the  whole  shaft  instead  of  a  single  ply,  three 
dimensional  analysis  of  plys  including  transverse  shear 
stresses,  delamination  problems,  thermal  effects,  and  others. 
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The  field  of  composites  also  opens  up  another  group  of 
considerations  like  the  effects  of  stacking  different  ply 
orientations  at  different  sequences,  unsymmetrical  laminates, 
unequal  thicknesses  of  plys  and  other  factors  relevant  only 
to  composites. 
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APPENDIX  A 


The  FORTRAN  program  used  in  this  investigation  is 
described  here.  Subroutine  ANALIZ  provides  the  basic 
analysis  used  in  the  optimization.  It  reads  the  initial 
design  description  and  calculates  the  values  of  the  objec¬ 
tive  function,  constraints,  and  all  other  parameters 
necessary  to  solve  the  analysis  problem.  COPES-CONMIN 
updates  the  design  to  minimize/maximize  the  objective 
function,  iterating  until  no  further  improvement  in  the 
objective  function  is  possible  without  violating  one  of 
the  constraints.  This  process  is  shown  in  the  flow  diagrams. 

The  global  catalog  lists  the  location,  FORTRAN  name, 
mathematical  symbol  and  description  of  the  parameters  used 
in  the  optimization.  These  parameters  are  contained  in  the 
labeled  COMMOM  block,  GLOBCM. 

There  are  at  least  five  input  cards  to  describe  the 
initial  design.  The  first  card  contains  the  initial 
dimensions — inside  radius,  thickness  and  shaft  length.  The 
second  card  describes  the  material  used.  Young's  modulus, 
shear  modulus,  Poisson's  ratio,  weight  per  unit  volume, 
yield  strength,  ultimate  strength  and  compressive  strength 
are  read  here.  The  endurance  limit  is  internally  calculated 
using  the  endurance  limit  factors  read  from  the  third  card. 
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The  size  factor,  CB,  is  internally  updated  as  the  shaft 
diameter  is  varied  during  optimization. 

The  fourth  and  fifth  cards  specify  the  design  and  load 
conditions.  The  parameter,  LC,  allows  the  user  to  choose 
the  end  conditions  applicable  to  his  requirements.  Pinned- 
pinned  and  clamped-clamped  analyses  are  available  while 
clamped-pinned  and  cantilever  end  conditions  are  yet  to  be 
developed.  The  parameter  LCC  allows  multiple  loading  design. 
One  can  thereby  design  a  shaft  that  is  capable  of  supporting 
LCC  different  load  combinations. 

The  magnitude  of  mass  imbalance  may  be  read  as  a  specified 
distance  from  the  center  of  the  shaft,  EO,  or  as  a  specified 
fraction,  EC,  of  the  shaft  radius.  The  safety  factor 
generally  used  in  the  analysis  is  also  specified  in  the  fourth 
card. 

The  fifth  card  specifies  the  horsepover,  shaft  speed  in 
RPM,  internal  pressure  and  axial  load  for  each  loading. 
Additional  cards  are  used  up  to  LCC  number  of  loadings.  The 
torque  is  calculated  internally  in  inch-pounds.  The  follow¬ 
ing  table  summarizes  the  required  information  for  subroutine 


ANALIZ. 


TABLE  IV 
•INPUT  CARDS 


CARD 

INFORMATION 

FORMAT 

1 

RI,  TH,  SC 

8F10.2 

2 

E,  G,  PR,  SPWT,  I/S 

SU,  YC 

8E10.2 

3 

CA,  CE,  CC,  CD,  CE, 

CF 

8F10.2 

4 

LC,  LCC,  EO,  EC,  SF 

215,  3] 

5 

HP,  RPM,  P,  F 

8F10.2 

Subroutine  BEND  calculates  the  bending  moment,  deflection 
and  shear  force  at  the  critical  cross-sections  of  the  beam. 
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SUBROUTINE  ANALIZ  (ICALC) 
FLOW  DIAGRAM 


ICALC  =  1 

READ  INITIAL  DIMENSIONS 
READ  MATERIAL  PROPERTIES 
READ  DESIGN  CONDITIONS 
PRINT  ALL  INPUT 

ICALC  =  2 

CALCULATE  GEOMETRIC  PARAMETERS 

CALCULATE  OBJECTIVE  FUNCTION 

CALL  SUBROUTINE  BEND 

.  CALCULATE  MAXIMUM  DEFLECTION 

.  CALCULATE  BENDING  MOMENTS 

.  CALCULATE  SHEAR  FORCE 

CALCULATE  STRESSES  AT  CRITICAL  STRESS 
ELEMENTS 

CALCULATE  CONSTRAINTS 

ICALC  =  3 

PRINT  ANALYSIS  RESULTS 


OPTIMIZATION  FLOW  DIAGRAM 
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TABLE  V 

GLOBAL  CATALOG  OF  PARAMETERS 


Global 

Location 

FORTRAN 

Name 

Math 

Symbol 

Definition 

1 

R1 

r . 

i 

Inside  Radius 

2 

TH 

t 

Thickness 

3 

SL 

i 

Shaft  Length 

4 

VOL 

V 

Volume 

5-14 

TQ(10) 

T 

Torque 

15-24 

P  (10) 

P 

Internal  Pressure 

25-34 

F  (10) 

F 

Axial  Load 

35-44 

RDS(10) 

w 

Shaft  Speed 

45 

E 

E 

Young ' s  Modulus 

46 
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Shear  Modulus 

47 

PR 
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Poisson's  Ratio 

48 

YS 

Sy 

Yield  Strength 
(Tension) 

49 

YC 

Sc 

Yield  Strength 
(Compression) 

50 

SU 

Su 

Ultimate  Strength 

51 

RO 

ro 

Outside  Radius 

52 

SPWT 

Y 

Weight  Per  Unit 
Volume 

53 

AI 

I 

Moment  of  Inertia 

54 

AJ 

J 

Polar  Moment  of 
Inertia 

55 

HS 

Hoop  Stress 

56 

AST 

°xA 

Axial  Stress  at  A 

57 

ASTN 

°x  (min) 

Minimum  Axial  Stress 

58-67 

ASC (10) 

°XB 

Axial  Stress  at  B 

68-77 

SALT (10) 

°a.  V 

Alternating  Axial 
Stress 

78 

SMEAN 

Mean  Axial  Stress 

79 

SST 

XT 

Torsional  Shear 
Stress 

66 


80 


SSV 


81 

82 

83 

84-93 

94-103 

104-113 

114 

115 

116 

117 

118 

119 

120-129 

130-239 


SSMX 

TCRT 


'max 


scsx 


VMS I (10) 

VMS 2 (10) 

VMS 3 (10) 

WM 

BM 

BMN 

V 

FC 


m 

°o- 

°mid 

CO 

c 

M 

max 
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SF  SF 
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Shear  Stress  due  to 
Sheer  Force 

Total  Shear  Stress 

Critical  Torque  for 
Buckling 

Critical  Stress  in 
Buckling 
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Von  Mises  at  x  =  0 
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APPENDIX  B 


The  derivation  of  the  equations  for  the  deflection, 
bending  moments  and  shear  force  at  the  critical  cross- 
sections  of  the  shaft  are  presented  here.  Pinned-pinned 
and  fixed-fixed  end  conditions  are  analyzed. 

1.  Pinned-Pinned 


_  -  m  HA 

Let  y  =  d  sm  -r- 

O  x. 

dQ  =  dM  (y+e) 

2 

dQ  =  pAw  (y+e)  dx 
Let  =  pAw^ 

dQ  =  Kx  [dQ  sin  +  e]  dx 


(1. 


(1. 


4 

75 


(1. 


With  the  effect  of  weight: 


dQ  =  [K1(dQ  sin  ^  +  e)  ~  K2]  dx 


(1.4) 


where  K2  =  pAg  =  numerically  equal  to  weight  per  unit  length 


-2  V 


Q'“  dQ  =  Kxr  -J-+eA]  -  K2  A 


(1.5) 


Moment  at  x  =  x 


M  =  -  *  (  Z  -  X  )  + 
xo  2  o 


I 


(x  -  XQ)  dQ 


(1.6) 


M  = 


.  Q  .  K 

2  +  K1 


9  d  _  d 

£  (~r  +  l}  +  [2  Ki£  (^r  +  e)]x  (1*7) 


Kie  T  0  -> 

+  4"x  -K1  do(t>  sinX 


K-  0  K, l 

-  T  x  +  4*  x 


This  represents  the  moment  at  any  value  of  x.  If  the  axial 
load,  F,  is  included,  the  total  bending  moment  is 

"TOTAL  *  Mx  +  **F  "  EI*" 


where 


Mp  -  Fy  =  F  (dQ  sin 


(1.8) 


Integrating  twice : 

K  © 

Ely  *  -i-  [/x  -  2 Ax3  +  x4]  +  Kx  dQ  (£) 4  sin  ~  d-9) 


+  F  dQ  (^)  2  sin  X  “  2?  t*3*  “  21ix^  +  ^ 
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At  x  =  l/0,  y  =  d 


o  4  o  2  3  K,  ei  5  Kj  Z 

Eld^  -  K.  d  (-)  -  Fd  (-)  =  — - - =4=- 

o  1  o  7T  0  7 r  384  384 


,  5£ 

ao  “  384 


K.e  -  K~ 

[ - ~  £  "4  - - TT^ 

El  -  K  (£)4  -  F(|)^ 


(1.10) 


This  is  the  deflection  when  the  weight  effect  is  subtractive. 
The  maximum  deflection  occurs  when  the  weight  tends  to  increase 
the  deflection. 


5£4  ^le  +  ^2 

max  384  ^EI  _  K  (A) 4  -  f(-) 

1  TT  IT 


(1.11) 


Substituting  (1.10)  in  (1.7)  and  (1.8),  the  bending  moment 


at  x  *  i/2  is 

Mmin  *  K1  1|+  do  <?»21  +  Fd  ' 


(1.12) 


Again  the  maximum  moment  occurs  when  the  weight  effect  is 


additive 


e  Z  2  ^2^2 

M  =  K,  [~  +  d  (-)  1  +  Fd  +  -4 — 

max  1  1 8  max  tt  j  8 


(1.13) 


From  Equation  (1.5),  the  maximum  shear  force  is  equal 

to  |  or 


K.  d^  Z  K.eZ  K-£ 

„  _  1  max  ,1  ,2 

VMAX  '  tt  2  2 


(1.14) 


2 .  Clamped-Clamped 

A  similar  derivation  procedure  for  clamped-clamped  ends 
using  a  shape  function: 


d 

O  n  2ttx. 

y  =  T  [1  -  cos  -J-] 


(2.1) 


and 


d  d 

dQ  =  k1 


~2~  cos  yy  +  e]  dx 


-£ 


Q  =  )  dQ  =  [-f  +  e] 

o 


(2.2) 

(2.3) 


Moment  at  x  =  x 


M  =  M  -  ^  x  + 
o  2 


f 


(XQ  -  x)  dQ 


K.  SL  d  (^q  x  d  x,d 

*  Mo  "  —  (T  +  e)  +  K1  j  ^  2 - 2  COS  T~ 


xd 

+  xoe  -  ~T 


xd 

~2 


o  2ttx  t  . 

cos  —x—  -  xe]  dx 


(2.4) 


M  at  any  value  of  x: 

K.  I  d  K  d  2 

Ely"  =  Mx  =  Mq - {-y  +  e)  x  +  -y  (-y  +  e)  x 


W  „  2ttx. 

~~2~  (1  -  cos  —» 

8tt 


(2.5) 
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T1 


Integrating  once: 


K1  do  x3  x2? 

Ely'  =  X  -  -f  (-£  +  e)  -  2^) 


K.  d  £  x 
1  o 

8^1 2 


K..d  £ 
1  o 


16tt 


y '  =  0  at  x  =  0 ,  2,/2  and  £,  ci  =  0 

Mo  -  K1  *2  [sf  +  nr  + 

8  TT 

Integrating  once  more: 


EiZ  =  fdo  ^  ^  ,x4  _  x3£  ,  xV 


K, 


+e) 


12 


24 


d  2. 

o 

32rr ' 


cos  +  C2 


y  = 


0  at  x  -  0,  C2  = 


15*  ,  <1° 

Ki 


d  £ 
o 

32tt' 

2 


(  2  +  e)  (I?)  (x  _  ^  +  ITT 


(l  - 


2  TTX\  /  ~  \ 

cos  —5—)  (2.6) 


32tt 


y  = 


dA  at  x  =  £/- : 

O  2 


__  d  „2  „2  d  £ 

El  ^  _  ,0  ,  ,£ \  ,  o 

Kx  do  (  2  e)  96  (  4  }  16it4 


i  i 


Eld  d  £4  .4  d  £4 

o  o  e£  o 

K,  768  384 


16tt 


4 
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o  = 


K1eS, ' 


KT  1 

384EI  -  -  (^y-  +  24) 


(2.7) 


With  axial  load  F: 


Fd 


Moment  due  to  F  =  M„  =  M  -  F  =M  -  — ~  (i  -  cos 

F  o  y_  o  2  i 

r 


Integrating  once: 


(2.8) 


EIyF'  =  mo  x  -  !f°  (x  -  i  sin  ♦  Cj 


y  '  =  0  at  x  =  0 ,  C,  =  0 

r  X 


Fdc 

yp  -0atx-e,Ho-T 


where  M  =  Moment  at  the  ends  due  to  the  axial  load 
o 


Fdo  ,x2 


Fd  2 

O  rX 


EIy  “  2^  (  2  ]  2  ^  2  +  4tt2  CoS  l 


2ttx 


]  +  C. 


y„  ~  0  at  x  =  0 ,  C-  = 


Fd  i 
o 

2  "  Q  2 

8tt 


Fdo£  „  2ttx^ 

y  =  - - -  (1  -  COS  -J-) 

F  8tt2EI  £ 


(2.9) 


At  X  —  Z/-  : 

Fdl2 

y  =  j 

4ttEI 
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From  Equation  (2.6) 


= 


K.  d  2 

1  r  O  ,X 


w  El 


[—  %>  (X  -  £)z  +  e  (^f)  (X  -  £) 


24' 


d  £' 

o 

32n^ 


(1  -  cos  2™)] 


At  x  =  Z/ 2 : 


y,.,  = 


K1^ 


w  384EI  1  2 


[-T  +  e]  + 


K.d  £ 

1  o 

16tt4EI 


d  =  y_  +  y 
o  1 F 


.  (j^L  +  V 


K1£ 


4  tt2EI  768EI  16,4EI 

2  K.£4  24K. £4 

384EI  d  =  +  _L_  +  - 1 

o  ,2  4 


d  = 
o 


2  TT 

4 

K1e£ 


d  + 


4 

K  £  e 
3  84  El 


)  d0+  K2  e£ 


384EI  -  Kx  £4  (|  +  -  96-F2£' 

IT11  IT2 


When  the  effect  of  weight  is  considered: 


do  = 


4  4 

Kxe£  ±  K2£4 

384EI  -  0.74638K1£4  -  96(|)F 


where  =  pAg,  numerically  equal  to  weight  per  unit 


(2.10) 


(2.11) 


(2.12) 


length. 
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